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ABSTRACT 

We consider the issue of selecting parameters and their associated ranges for carrying 
out searches for continuous gravitational waves from steadily rotating neutron stars. 
We consider three different cases (i) the ‘classic’ case of a star spinning about a prin¬ 
cipal axis; (ii) a biaxial star, not spinning about a principal axis; (iii) a triaxial star 
spinning steady, but not about a principal axis (as described in Jones, MNRAS 402 , 
2503 (2010)). The first of these emits only at one frequency; the other two at a pair 
of harmonically related frequencies. We show that in all three cases, when written in 
terms of the original ‘source parameters’, there exist a number of discrete degeneracies, 
with different parameter values giving rise to the same gravitational wave signal. We 
show how these can be removed by suitably restricting the source parameter ranges. In 
the case of the model as written down by Jones, there is also a continuous degeneracy. 
We show how to remove this through a suitable rewriting in terms of ‘waveform param¬ 
eters’, chosen so as to make the specialisations to the other stellar models particularly 
simple. We briefly consider the (non-trivial) relation between the assignment of prior 
probabilities on one set of parameters verses the other. The results of this paper will 
be of use when designing strategies for carrying out searches for such multi-harmonic 
gravitational wave signals, and when performing parameter estimation in the event of 
a detection. 

Key words: gravitational waves - methods: data analysis - stars: neutron - stars: 
rotation 


1 OVERVIEW 


Rotating neutron stars are potentially detectable sources of 
continuous gravitational radiation. They may emit a steady 
gravitational wave signal because of a non-axisymmetry in 
their mass distributions, caused either by elastic strains in 
their solid phase(s), or by magnetic strains sourced by their 


global magnetic field (see e.g. Andersson et al. (20111 for a 


review). A star with such a ‘mountain’ will typically emit 
at a frequency 2/, where / is the spin frequency. If the star 
does not rotate steadily but instead undergoes free preces¬ 
sion, the gravitational wave signal is then emitted at fre¬ 


quencies equal (or close to) both / and 2/ (Zimmermann 
fc Szedenits|1979 Jones fc Andersson|2002 I. However, such 
precession would normally leave an imprint on the observed 


radio pulsation received from a pulsar (see e.g. Jones & An 
dersson (20011). Such modulations are not typically seen 
in the known pulsars. For this reason, most gravitational 
wave searches to date that have targeted known pulsars have 


searched only at the frequency 2/; see Aasi et al. (2014 1 and 
references therein. (The exceptions have been two ‘narrow 
band’ searches for the Crab pulsar and one for the Vela pul¬ 


sar, where a small band around 2/ was searched (Abbott 
|et al.| ( [2008[ ), and |Aasi et al.| ( |2015[ )). 


However, as shown in Jones (20101, the presence of a 


pinned superfluid within the star can change this picture. A 
star can then rotate steadily, and still produce gravitational 
radiation at both / and 2/, providing the axis about which 
the pinning takes place does not coincide with a principal 
axis of the star’s moment of inertia tensor. This motivates 
the carrying out of searches for such multi-harmonic signals 
from known pulsars, despite their lack of precession. 

Given these considerations, we can identify three types 
of continuous gravitational wave emission. There is the gen¬ 


eral case (as per Jones (2010l), which we term the triaxial 


non-aligned case, with emission at both / and 2/. There 
is also the simplest case, of a triaxial star spinning about 
a principal axis. This is the sort of signal assumed in most 
targeted gravitational wave searches to-date, and produces 
emission at only 2/. We term this the triaxial aligned case. 
There is also an intermediate biaxial case, where the star 
is assumed biaxial. This produces gravitational radiation at 
both / and 2/, with a waveform of intermediate complexity. 
The waveform in this case is identical to that of a biaxial 
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precessing star, as considered in the literature (Zimmermaiin 
fc Szedenits|1979j Jones fc Andersson|2002 1 , SO is an impor- 


tant case to include, but note that in the model of [Jones] 
(20101, the gravitational wave emission is not accompanied 
by precession. 

A number of parameters appear in these models; we 
term these conventional parameter choices the source pa¬ 
rameters. It turns out that there are two issues with these 
parameters, which we address in this paper. Both issues re¬ 
late to the existence of degeneracies, i.e. the existence of 
different values of the source parameters that produce the 
same detected gravitational wave signal h{t). Firstly, for all 
of the models, there exist a number of discrete degeneracies. 
These are related to discrete symmetries of the star’s mass 
distribution. Secondly, in the case of the model and param- 
eterisation of Jones (20101, there also exists a continuous 
symmetry, i.e. there exists a 1-parameter family of source 
parameters that all produce the same h{t). 

We do two main things in this paper. We first identify 
the discrete symmetries of the models written in terms of 
source parameters. This allows us to give minimal ranges 
that are needed in these parameters, removing redundancy 
in this form of the parameterisation, such that there is a 
unique set of parameters corresponding to any given star. 
We present such ranges in tabular form. Secondly, we iden¬ 
tify a convenient set of waveform parameters, in which the 
continuous degeneracy present in the parameterisation of 
|Jones| ( |2010[ ) is removed by effectively reducing the number 
of parameters by one. Any future gravitational wave search 
could be first carried out using this waveform parameter¬ 
isation. The conversion to the corresponding 1-parameter 
(but possibly more insightful) family of source parameters 
could then be carried out by making use of formulae given 
here, with the ranges in source parameters restricted appro¬ 
priately. We also (very briefly) discuss the relationship be¬ 
tween prior probabilities assigned to the source parameters 
and the corresponding prior probabilities on the waveform 
parameters. 

The existence of the continuous degeneracy of the model 
of Jones (20101 was noted by Bejger & Krolak (20141, who 
wrote down a reduced parameter set that removed this de¬ 
generacy. The waveform parameters used here are differ¬ 
ent from the parameter set introduced by [Bejger fc Kr61ak| 
(20141, and more closely tied to the fundamental scalar 
quantity, the mass quadrupole moment, that described the 
gravitational wave emission properties of the star. Our cho¬ 
sen parameterisation makes the specialisation to simpler 
forms of gravitational wave emission particularly transpar¬ 
ent. 


The results of this paper will be of use to gravitational 
wave observers when devising strategies for carrying out 
searches for gravitational wave signals with such multiple 
frequency components. A study of the issues raised by such 
searches is currently underway, and will be presented else¬ 
where (Pitkin et ah, in preparation). 

The structure of this paper is as follows. In Section 
we look at the source parameter description, carefully iden¬ 
tifying the discrete symmetries, and thereby finding minimal 
ranges to which the parameters can be restricted. In Section 
[^we give the waveform parameterisation, in which the con¬ 
tinuous degeneracy is removed. In Section]^ we briefly dis¬ 
cuss the relationship between priors written in terms of the 


source parameters and priors written in terms of the wave¬ 
form parameters. We summarise our findings in Section 
For the convenience of those carrying out gravitational wave 
searches, in the Appendix we provide summary tables show¬ 
ing (one possible choice) of ranges appropriate to the source 
parameters (Table [CT| and the waveform parameters (Table 


2 CHOOSING RANGES IN THE SOURCE 
PARAMETERS 

We will consider gravitational wave emission from a star 
spinning steadily, i.e. with an angular velocity ft, fixed in the 
inertial frame. The star has a moment of inertia tensor that 
is constant in the rotating frame, with principal components 
{Ii, I 2 , Is). In the most general case, the rotation axis need 
not coincide with one of these principal axes; as argued in 


Jones (20101, the presence of an internal pinned superfluid. 


with pinning axis misaligned with a principal axes, will be 
of this class. As described above, and as will be elaborated 
upon below, there are two special cases, the triaxial aligned 
case and the biaxial case. There are a number of features 
common to all three cases, which we will now describe. 

The signal h{t) received by a gravitational wave detec¬ 
tor is given by equation (4) of jjaranowski et al.j ( |1998[ ) , here¬ 
after JKS: 


h{t) = F+{t)h+{t) + F^{t)h^{t), (1) 

where the antenna functions F+, Fx are given by equations 
(10) and (11) of JKS: 

F+{t) = sinC[a(t) cos2i/)poi-I-&(t) sin2V>poi], (2) 

Fx{t) = sm([b{t)cos2ippo\ — a{t)sin2ippoi], (3) 

where ^ is the angle between the interferometer arms and 
ippoi is the polarisation angle; the subscript ‘pol’ has been 
added to avoid confusion with the Euler angle tp that ap¬ 
pears below. The functions aft) and b{t) are complicated 
functions of source location in the sky (specified by right 
accession a and declination <5) and detector location on the 
Earth, as given by equations (12) and (13) of JKS. The 
phasing of the signal, as give by JKS equation (14), also 
depends upon the source location in a complicated way, so 
we will always need to cover the full sky parameter space of 
0 < a < 27r, —7r/2 < <5 < 7r/2. Physically, the spin vector 
can point in any direction, over the full range in inclination 
angle 0 < t < tt and the full range in polarisation angle 
0 < V'pol < 271. 

Collecting these results, we see that in all cases, the 
following ranges in gravitation wave frequency, sky location 
and spin vector orientation correspond to physically distinct 


sources: 




0 < 

Q, 

< 00, 

(4) 

0 < 

a 

< 2tt, 

(5) 

-7r/2 < 

6 

< 7r/2, 

(6) 

0 < 

i 

< tt, 

(7) 

0 < 

V^pol 

< 277. 

(8) 


Strictly, n is a function of time, so = flft). In practice, 
the full phase evolution is often parameterised as a Taylor 
expansion in fl and its time derivatives; by writing the single 
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2: 



Figure 1. The orientation of our body is specified by the three 
standard Euler angles (0, </>,-i/;), as labelled above. The fixed 
inertial-frame axes are denoted by (x, y, z), while the body-frame 
axes are denoted by {x,y,z) and rotate about the inertial 2 :-axis. 
The so-called line of nodes, N, lies along the intersection of the 
xy and xy planes. 


parameter Q, we are subsuming all such phase information 
into the one parameter, to avoid introducing further param¬ 
eters that have no bearing on our considerations here. 

Given the form of equations § and if one only 
cares about the received waveform h{t), it is clearly possible 
to restrict the range in tppoi further: 

0 < tppol < TT. (9) 


In fact, as we argue below, if one only cares about the re¬ 
ceived waveform h{t), and has no ‘prior’ information on the 
other parameters, it is possible to restrict the range in ippoi 
further to 0 < i/ipoi < n/2. However, the full range over a 
complete circle given in equation Q is the one required to 
represent all physically distinct stellar configurations, and 
a value for tppoi obtained by non-gravitational wave means 
could lie anywhere in this interval. 

We now turn to the form of the polarisation components 
h+{t) and hx{t). In the model of Jones (20101, the orienta¬ 
tion of the body, and therefore the gravitational wave emis¬ 
sion, depend upon three angles 9, 4>, tp, basically just Euler 
angles giving the orientation of the body with respect to 
the inertial frame; see Figure (This triaxial non-aligned 
case has the triaxial aligned and biaxial solutions as special 
cases). In the pinned superfluid case 9 and ip are constants, 
while (p Is the angle that generates the rotation, so that 


(p — lit (pQ , 


( 10 ) 


with (po a constant, giving the orientation of the body at 
time t = 0. However, as we show below, the phase function 
that actually appears In the waveforms is c^gw (or twice this), 
given by 


= Qt + ( 


^gw,0? 


( 11 ) 


where we have defined the constant 


0 gw,O — 00 0obs5 (12) 

where 0obs is the azimuthal location of the observer. It was 
obvious that a constant of this form should appear In the 
waveform, as it is only the t = 0 position of the source 
relative to the observer that can affect the received signal. 

Also, it is only the asymmetries in the moment of inertia 
tensor that appear in the waveforms, so we define: 

A721 = I2 — I1, (13) 

AI31 = I3 — Ii. (14) 

We therefore see that we have a set of ten parameters, 
which we refer to as the source parameters: 

^source — {^5 ^1 V^poli ^-^ 217 ^-^31 , 0gW,O5 Ip}. (15) 

The ranges in the hrst hve parameters that correspond to 
physically distinct stellar configurations were given in equa¬ 
tions above. 

The ranges and values of the other parameters depends 
upon which of the three cases we consider, but their maximal 
ranges are: 

0 < 0 < TT, 0 < 0gw,o < Stt, 0 < Ip < 2' k , (16) 

for the Euler-type angles. For the sake of definiteness, we will 
also assume there is some restriction on the allowed sizes of 
the asymmetries in the moment of inertia tensor, although 
this really depends upon the (poorly constrained) physical 
mechanisms that produce the deformation in the first place: 

— A/max < A/21 < A/max, —A/niax < A/31 < A/max. 

(17) 

Some of these parameters can be set to zero, or simply don’t 
appear, in the triaxial aligned and biaxial waveforms, while 
the appropriate ranges are also dependent upon the model. 

In this Section, our purpose is two-fold: 

(i) To identify particular ranges in the parameters 
{A/ 21 , A/ 31 , 0, 0gw,o, ”0} such that all physically distinct 
configurations of a star’s mass quadrupole (which is the 
quantity responsible for generating the gravitational waves) 
can be described uniquely. It will be possible to incorporate 
all conceivable additional prior information on the param¬ 
eters (possibly obtained by electromagnetic means) within 
these ranges. This will eliminate most of the degeneracies in 
the waveform. We provide a summary of a possible choice 
of source parameter ranges in Table m 

(ii) To identify and exploit one further discrete symme¬ 
try in the waveform connected with the polarisation angle 
0poi that, in the absence of prior information, would allow 
a search to be carried out over a slightly smaller parameter 
range (leading to a simpler search), together with a rule for 
generating the other, equally acceptable parameter values. 
In the event of a successful detection. 


In terms of our first aim, we will exploit the fact that, 
as we are considering gravitational waves generated by the 
mass quadruple of the star, we only care about the actual 
orientation of the star up to a tt rotation about any one of 
the body axes (x, y, z), rendering the full ranges of equation 
(161 redundantly large. For instance, suppose some glowing 
hotspot is observed in the rotational equator of a triaxial 
aligned star, and some astronomer’s theory said that this 
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must correspond to the position of the axis of least moment 
of inertia. Suppose this hotspot is such that at time t = 0 
it lies on the far side of the star, relative to the observer. 
We would then set <?!>gw,o = tt. However, the tt rotation sym¬ 
metry of the mass quadrupole means that there must also 
be an axis of least moment of inertia on the nearside of the 
star at t = 0, so we could equally well set (?i>gw,o = 0, i.e. 
we have the freedom to map all values of </)gw,o into the 
range 0 < (f>gw,o < tt- Other rotational symmetries can be 
exploited for the biaxial and triaxial non-aligned cases, al¬ 
though their description in terms of Euler angles will be 
more complicated, as will be described below. 

In terms of our second aim, note that it is obvious from 
the equations of JKS (i.e. equations above), a ro¬ 

tation in orientation angle 'i/ipoi —>■ 'i/’poi + ^/2 produces a 
change in sign of h{t), i.e. h{t) —^ —h(t). But there is a sec¬ 
ond way of producing an (also physically distinct) star with 
waveform —h(t). Consider emission from a star described by 
a given set of parameters, producing a waveform h{t). Now 
consider emission for a star that has an identical spin vector 
and sky location, but whose density perturbation 5p away 
from sphericity is reversed in sign, i.e. 5p —> —5p. Such a 
star will produce a waveform —h{t), and will be described 
by a different set of parameters. If follows that if both op¬ 
erations are carried out at once the waveform is unchanged, 
i.e. there is a degeneracy in the waveform. Explicitly, given 
a signal with a particular set of parameters, there will exist 
three other sets of parameters corresponding to the same 
signal, obtained, by transforming tppoi —>■ V’poi + 7 r /2 and 
simultaneously transforming some or all of the parameters 
{ 0 , <(>gw,o, V’, A/ 31 , A 72 i}, in a way that depends upon the 
choices already made for the allowed ranges of these param¬ 
eters. Note that solutions that differ by tppai —>■ tppoi + tt/2 
will be described by different { 0 , ())gw,o,'i/’i A/31, A 72 i} pa¬ 
rameters, while solutions that differ by tppoi —>■ ippoi + "'iU 
be described by the same values of { 6 , <)>gw,o, V’, A/31, A/21}, 
so differ only in the orientation of their spin vectors. The 
existence and effects of this sort degeneracy were described 
long ago in the context of biaxial processing stars by |Zim-| 


mermann & Szedenits (19791. 


This means that there several ways in which the polar¬ 
isation angle can be handled in a gravitational wave search: 


(i) If electromagnetic observations have provided a value 
(or at least small range) in i/jpoi (whose value may lie any¬ 
where in the range 0 < V'poi < Stt), then this should be used 
in the search, and a detection will correspond to a single set 
of parameters {A/21, A/31, 0, <()gw,o, V’}- 

(ii) If electromagnetic observations do not constrain ippoi, 
one could search over the interval 0 < ^poi < 27r, but 
this would be highly redundant, with a detection potentially 
manifesting itself at four different points over the searched 
parameter range, with each inferred parameter set differing 
by tv/2 in i/'poi. 

(iii) More sensibly, one could search over the reduced 
range 0 < ippoi < tt, but again mindful of a degeneracy, 
with a detection manifesting itself at two different values of 
Tppoi within this range, differing by tt/2, and corresponding 
to two different parameter sets {A/ 21 , A/ 31 , 6, (pgw,o, '!/’}■ The 
remaining parameter sets, corresponding to identical wave¬ 
forms, can then be generated by the simple transformation 
Tppoi —>■ V'pol + TT, keeping the other parameters fixed. 


(iv) Probably most sensibly of all, one could search over 
the reduced range 0 < Tppoi < 7 r/ 2 . This eliminates the 
discrete degeneracy, and a detection would be expected to 
manifest itself as a single set of parameters. In the event of 
a detection, three other sets of parameters, corresponding 
to identical waveforms, can then be generated by successive 
use of the transform i^ipoi —>■ i/’poi+ 7 r/ 2 , while simultaneously 
carrying out a transform on (all or some) of the parameter 
set {A/21, A/31, 0, ((>gw,o, V’li in 8 - "'ny that depends upon 
the model under consideration, as will be described in this 
Section. 


We will now turn to a consideration of each of our three 
models. When we write out the polarisation components h+ 
and /lx below, we refer to a gravitational wave detector with 
its 1 -arm along eg and its 2 -arm along relative to the in¬ 
ertial frame axes (x,y,z) described above. See Jones (20121 
for the generalisation to an arbitrary detector orientation. 


2.1 Triaxial star, not spinning about a principal 


In the general case the wave field can be written as (see 
Jones (20101 and Jones] ( [2012[ ) for details): 


y2n _ ^^||^/2i(sin^-i)! — cos^ i/) cos^ 0) 

—A /31 sin^ 6] cos 2<)>gw + A /21 sin 2Tp cos 6 sin 

20 212 ^ r 

/lx =-2 cos i{ A /21 sin 2'(/) cos 6^ cos 2()igw 

r 

— [A /21 (sin^ Tp — cos^ Tp cos^ 0) 

—A /31 sin^ 6] sin 2(pgvt}, 


T 

/l I = 


rll 

= 


sin i cos i{ A /21 sin 2Tp sin d cos (p. 


gw 

-|-(A/ 2 i cos^ Tp — A/ 31 ) sin20sin())gw}, 
r 


— — sinij (A /21 cos^ Tp — A/ 31 ) sin20cos<()g 
r ^ 

—A /21 sin 2tP sin 0 sin <))gw }. 


(19) 


( 20 ) 


( 21 ) 


(This differs only trivially from the waveform given in|Jones 


( 20101 , were the observer location was fixed to (poha = —tt 12, 
the detector with respect to which and hx were referred 
had its 1 -arm along = eJ, and its 2 -arm along = —e®. 


See Jones (20121 for the equations for the general metric 


perturbation h^b, and for projecting this onto an arbitrarily 
orientated detector). 

The physical and therefore default ranges in the Euler 
angles are given by equation ( |16[ ) above. One could allow 
the parameters A /21 and A /31 to take either sign (positive 
or negative). However, we are free to follow the common 
convention of rigid body dynamics and choose our axes such 
that 13 > I 2 > /i, so that 


A/,, 


> A /31 > A /21 > 0. 


( 22 ) 


There exist discrete symmetries that can be exploited 
to reduce the range in parameters further. The quadrupole 
moment tensor, and therefore the gravitational wave field, 
is invariant under rotation of tt about any one of the three 
body axes. Ox, Oy or Oz . A rotation about Oz corresponds 
to Tp ^ Tp + TT (this is obvious from the definition of Tp, 
but a proof is given in Appendix A1 1 . The waveform above 
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clearly only depends upon sin 2 'tp and cos 2 ip (or as functions 
that we be re-written in terms of these), so indeed has this 
symmetry. This means we can halve the range in this angle 
to 0 < V' < ^, leaving 

0 < 6 » < TT, 0 < <()gw,o < Stt, 0<1p<TT. (23) 

A rotation of tt about the body’s Oy axis corresponds 
to the mapping ( 6 , ())gw,o, V’) —>■ (tt — '*/’); see 

Appendix I A2[ The waveform above can indeed be shown to 
posses this symmetry. This means we can make a further 
reduction, halving the parameter range of any one (but only 
one) of the parameters (9, ((>gw,o, V’)- (Note that the waveform 
must also be invariant under a tt rotation about Ox, but this 
is equivalent to the composition of the above two rotations, 
so cannot generate any further reduction in the parameter 
space.) This means there are three options for the Euler 
angle parameters: 


0 < 9 < 7r/2, 0 < 4>g-w,o < Stt, 0 < ip < tt, or (24) 
0 < 9 < n, 0 < <))gw,o < tt, 0 < V’ < tt, or (25) 
0 <9<n, 0 < (?igw,o <271, 0 < ip < 7r/2. (26) 

Any one of these choices, together with the ranges of equa¬ 
tion (221 and the parameter ranges of equations will 

always be able to accommodate any additional prior infor¬ 
mation on the star. 

However, this parameter space is redundantly large 
from the point of view of carrying out a gravitational 
wave search without additional prior information. As ar¬ 
gued above, if there exists a star producing a wavefield with 
components there must exist another, physically dis¬ 

tinct star, whose wavefield has the sign of all these compo¬ 
nents reversed, i.e. Physically, this cor¬ 

responds to reversing the sign of the density perturbation 
5p that deforms the star away from spherical symmetry. 
The transformation that produces the mapping —>■ 

is rather complicated. The details are given in Ap- 
pendix[^ and involve a transformation mixing the Euler an¬ 
gles {9,(pg-„fi,ip) and also mixing the amplitude parameters 
(A721, A/31). If this transformation is carried out together 
with the operation ^poi —>■ ippo\ + 7r/2, the full waveform 
h(t) is invariant. The significance of this is that a further 
reduction in the parameter space is possible. One option is 
to reduce the range in ippoi to the range: 


0 < V’poi < 7r/2. 


(27) 


There is presumably another option, involving some reduc¬ 
tion in the parameter set {9, <j)g,„fi,ip, AI 21 , AI 31 ), but it is 
not clear how to implement this second option. For instance, 
the Euler angles aren’t simply increased by tt or multiplied 
by — 1 , so the reduction presumably isn’t a simple halving 
their ranges (see equations ( B13[ |-( B151 in Section 0 . 

If a gravitational wave search is carried out over this 
restricted range 0 < ippoi < 7 r/ 2 , three other equally accept¬ 
able solutions can be obtained by successive applications of 
the transformation 


ippoi ->■ Ippol + 7r/2, 


(28) 


with 


corresponding transformation 


of 


(0, (/)gw,o,A/ 21 , A/ 31 ) given by equations (B13|-(B15l. (If 
one wishes the parameters to remain conhned to one of the 
minimal ranges identified above, further transformations 


may be needed, e.g. the tt rotation about Oz of Appendix 
All. 


Note that if we weren’t enforcing the inequalities of 
equation \ 22 \ and were instead allowing the quantities A /21 
and A /31 to take either sign, the operation Sp —^ —5p could 
be achieved much more simply, by making the replacements 
A /21 —>■ —A /21 and A /31 —>■ —A/ 31 . This would, however, 
lead to different and less easily derivable choices in minimal 
ranges of the Euler-type angles. 

See the first column of Table |C1| for a summary of a 
possible choice of parameter ranges, where, for definiteness, 
we use the ranges of equation (241 for the Euler-like angles. 


2.2 Intermediate case: A biaxial star, not 
spinning about a principal axis 


The simplest (and conventional way) of describing a biaxial 
body is to single ont the 5-axis as special, i.e. to set Ii = I 2 , 
so that A /21 = 0. The waveform can then be shown to be 
( Jones|2012 I 


i 2 Q, 

h+ = 

— ^^(1 + cos^ t)A/3i sin^ 6 cos 20gw, 
r 

( 29 ) 

, 20 . 

/lx 

^ 2 

- 2 cos tA/31 sin ^sin20gw, 

r 

( 30 ) 

7 r2 

h+ = 

- sin L cos lAIsi sin 20 sin 0 gw, 

r 

( 31 ) 

1 Q 

/lx = 

— sin lAIsi sin 20 cos 0 gw 

r 

( 32 ) 


This is the wave field of a biaxial star spinning about an axis 
other than a principal axis. In the model of Jones (20101 it 
corresponds to a non-precessing star with a pinned super¬ 
fluid. It is also identical to the GW field of a processing 
biaxial star without pinning, of the sort considered by |Zim-| 
[mermann fc Szedenits|p979p and JKS. So, it is an important 
case to cover. 

By singling out the 5-axis as the symmetry axis, we have 
A /21 = 0. Having made this choice, we can no longer insist 
that I 3 is the axis of greatest moment of inertia. Instead, 
we must allow for the star being either oblate (A /31 > 0) or 
prolate (A /31 < 0), depending upon the sign of A/ 31 , so we 
have: 


A/21 = 0, —A/max < A/31 < Alir 


(33) 


Now consider the Euler angles, whose ‘default’ ranges 


were given in equation (I 61 . We can again exploit symme¬ 


tries. The waveform no longer depends upon the angle ip, 
a consequence of the axisymmetry of the body about the 
Oz body axis, so this angle is removed from our considera¬ 
tions. As described in Section [A2| the operation of perform¬ 
ing a rotation of tt rotation about the Oy axis takes the form 
{9, <pgvi,o) —^ (tt-0, <pgvi,o + TT), allowing us to halve the range 
in 9 or <(>gw,o- We can therefore have: 


0 < 6^ < tt/2, 0 < <pgw,o < 271, or 

0 < 6^ < 71 , 0 < <(>gw,0 < TT. 


(34) 

(35) 


Either one of these two options, together with the parameter 
ranges of equation (331 and equations will always 

be able to accommodate any additional prior information on 
the star. 

However, this parameter space is again redundantly 
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large from the point of view of carrying out a gravitational 
wave search without addition prior information. To see this, 
note that the transformation A /31 —>■ —A /31 (equivalent to 
the transformation 5p —>■ —5p) changes the sign of the above 
polarisation components, i.e. thereby flip¬ 

ping the sign of h{t). The transformation ipp^i —>■ ippoi + 11/2 
also flips the sign of h{t), so the two transformations to¬ 
gether leave h{t) unchanged. It follows we can reduce the 
range in A /31 or the range in "i/jpoi, i.e. we have the choice 


0 ^ 'ippol ^ ^/2, /max ^ A /31 -^max ) or (36) 
0 < V'pol < TT, 0 < A /31 < -^max ■ (37) 


This degeneracy was noted by [Zimmermann fc Szedenitsj 
(19791. A gravitational wave search for a biaxial star could 
then be carried out with either of the two choices hardwired 
in, with the understanding that in the event of a detection, 
three other equally valid solution can be obtained via suc¬ 
cessive uses of the transformation: 


^poi ->■ tppol + 7 r/ 2 , A /31 ->• -A /3 


(38) 


See the middle column of Table |C1| for a summary of 
a possible choice in parameter ranges, where for definite¬ 
ness we choose the ranges of equation (341 for the Euler-like 
angles. 


2.3 Simplest case: A triaxial star, spinning about 
a principal axis 


This is the standard case of a triaxial star spinning about a 
principal axis, emitting only at 20 ., i.e. the sort of emission 
normally assumed in continuous gravitational wave searches. 
The convention is to choose the rotation axis to be the z- 
axis. This is accomplished by setting 6 ^ = 0 in equations 
(Tspri to give 


hjp = -A/2 i(1 -|-cos^ t) cos 2[flt-I-(())gw,o + V’)]>(39) 

r 

20 ^ 

hy = -A/2i2cosrsin2[flt -I-((()gw,o + V’)]- (40) 

r 


Note that the parameters angles (()gw,o and ip are degenerate, 
i.e. only their sum appears in the waveform. 

We are free to lay down the {Ox, Oy) axes such that 
A /21 > 0. The symmetry of rotating by tt about Oz then 
corresponds to {pgwfi + ip) — >■ (<(>gw,o +'*/’)+ '^- The waveform 
clearly has this symmetry, suggesting we need to cover the 
range 0 < (<(>gw,o+ '*/’) < As we have fixed 9 = Q, this rota¬ 
tion about Oz is the only angular degree of freedom, so there 
are no further symmetries we can exploit, corresponding to 
rotations about Ox or Oy. So, for a triaxial aligned rotator, 
the set of physically distinct configurations is spanned by: 


A /21 >0, 0 < (((>gw,o -I- V') < ^- (41) 


These choices, together with the parameter ranges of equa¬ 
tions will always be able to accommodate any ad¬ 

ditional prior information on the star. 

However, these parameter ranges are redundantly large, 
in the sense that, in the absence of such prior information, 
smaller parameter ranges can be used to carry out gravi¬ 
tational wave searches. To see this, note that the waveform 
changes sign {h{t) —>■ —h{t)) when either one of the following 
transformations is performed: 


(i) ippoi ^ ippoi T 7r/2 , 

(ii) (<()gw.o + V’) (<(>gw.o + ’P)+ 7r/2 • 

The second of these transformation is equivalent to swapping 
over the axes of largest and smallest moment of inertia that 
lie in the rotational equatorial plane; as such it is not the 
same as the transformation Sp —>■ —5p discussed above, but 
rather flips the sign of h{t) in a way that preserves our choice 
of fixing the Oz axis as the rotation axis. 

It follows that, in carrying out a gravitational wave 
search, we can reduce the range in any one (but only one) 
of these parameters, so the options are: 


0 < V’poi < 7r/2, 0 < {<pgw,o + V*) < TTi or (42) 

0 < V'poi < TT, 0 < {<pgw,o + 'ip) < 7r/2. (43) 


The first of these three options has traditionally been used 
in gravitational wave searches (see e.g. Aasi et al. (2014|), 
although it should be noted that the phase angle that ap¬ 
pears in the literature is actually 


4 >gw,o = 2 ((()gw,o + Ip), (44) 


so that searches have traditionally searched over the range 
0 < $GW,o < 27r. 

If one of these restricted parameter spaces is used and a 
signal detected with parameters (V'poi, {(pgvr.o+ip)), three ad¬ 
ditional solutions can be obtained by successive applications 
of the transform 


ippoi —>■ 'i/'pol+ 7r/2, ((^gw.o + V') {(f>gv/,o~\~ip)~\~Ti/2. (45) 

Without additional (non-gravitational wave information) all 
such solutions are equally valid. (If one wishes the parame¬ 
ters to remain confined to, say, the range 0 < {pgw.o + V') < 
TT, then a further transformation {<j)gw,Q+ip) —t ('(>gw,o+V’)+'^ 
can be applied when necessary). 

See the final column of Table [UT] for a summary of these 
possible choices in parameter ranges. 


3 REFORMULATING IN TERMS OF 
WAVEFORM PARAMETERS 

The 10-parameter triaxial non-aligned waveform of Section 
|2.1| contains a degeneracy, and in fact only depends upon 9 
parameters. An easy way of seeing this is to note that if the 
cosine and sine terms in each of the four equations giving 
the polarisation components (equations ( |18[ )- are com¬ 
bined into single trigonometric terms (essentially writing the 
equations in ‘amplitude-phase’ form), the five parameters 
(0, ((igw,o, V"! A/ 21 , A/ 31 ) appear in only four different com¬ 
binations. Another, possibly more insightful, way of under¬ 
standing this is to return to first principles, making use of 
the multipole formalism for gravitational wave emission, as 


described in Thorne (19801. 


The fundamental quantities that appears in the wave 
generation equations Thorne are the mass quadrupole mo¬ 
ment scalars, related to the source’s density field p by equa¬ 
tion (5.27a) of Thoriie] ( |1980[ ): 


r 2 m 167 r -\/3 / , 2 


jZTn _ 


15 


/ 


pYFr^rUV. 


(46) 


The transverse traceless (TT) description of the GW field is 
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given by equation (4.3) of |Thorn^ (19801: 


1 TT /, \ _ 1 \ ' T' 2 mr 7 -iE 2 , 2 m 

n^ab W - - ^ab 

m 


(47) 


where is a tensor spherical harmonic. 


For rigid rotation about the z-axes at rate Q, the mass 


quadrupole scalars can be shown to take the from (Jones 


20121 


j2m _ ^complex —im(r2t+^o) 

1 — Oo™ e , 


(48) 


where ^ complex number that encodes details of 

the source, and 4 >o is as defined above, i.e. a phase angle 
giving the rotational phase of the body at time t = 0. We 
can write in amplitude-angle form: 


where C 2 m = > 0 and 0 < $ 2 m < Ztt. 

The waveform for an arbitrary rigidly rotating source 
is then given by equation \A7\ . In writing it down, it is con¬ 
venient to include some additional factors in our amplitude 
parameters; we define 


<52 


m — 



(50) 


so that the waveform can then be shown to take the very 
simple form ( [Jones |[M^ : 


/i+(2H) = 

— C 22 COs[2nt -1- $^2](1 + COS^ t). 

(51) 

h^{2n) = 

—C 22 sin[2Ht -|- $^ 2)2 COS t. 

(52) 

/i+(H) = 

— ^(521 cos[nt -|- $^i] sintcost. 

(53) 

h^{Q.) = 

-^(521 sin[nt -(- sint, 

(54) 


where the phases ^ 2 m are related to previously introduced 
quantities by 

= 2 <^gw,o-$ 22 , (55) 

= <(>gw,o - $21. (56) 


These equations can then be specialised to the three cases 
considered above. They are clearly rather simple in form, 
with all the (potentially) complicated details of the source 
parameters being buried within the amplitudes C 22 and C 2 \, 
and the phases ^g and $g. 

This approach also has the advantage of making the 
counting of the number of parameters more straightforward. 
We can count as follows. We need the set of hve parame¬ 
ters {n, a, S, b, 'ippoij giving the spin frequency, sky location, 
and spin orientation of the source, as before. (The maximal 
ranges in these parameters were given in equations 
earlier). For a steadily rotating source emitting gravitational 
waves only at 2 Q, we then have the amplitude-phase pair 
C 22 ,^g also, giving seven parameters, consistent with the 
number of source parameters in this case. However, in the 
triaxial non-aligned case, where the H-harmonic is present 
too, we also have the amplitude-phase pair C 21 , $21 ■ This 
gives a total of nine parameters, not the ten that one would 
arrive at by examining the waveform as written previously, 
confirming the existence of a continuous degeneracy in the 
source parameters in this triaxial non-aligned case. For the 
biaxial case, we will find that there is a particular relation 
between the phases $21 and $221 giving eight parameters. 


equal to the number of source parameters, so there is no de¬ 
generacy for biaxial stars, only for triaxial non-aligned ones. 

We can collect the relevant parameters together to give 
the nine waveform parameters: 

^waveform — ^5 <■) ^pol, <521, $21 ; <^22 , $22 } • (^'^) 

Comparing with the ten source parameters of equation ( |15| ) , 
we see that the Hrst five parameters {fl, a, J, t, V'poi} are com¬ 
mon between the two parameterisations, while the set of five 
source parameters {A72i, A/ 31 ,0, (()gw,o, V’} are replaced by 
the set of four waveform parameters {< 52 i, $ 21 , C 22 , ^g}- 

It would therefore seem that there may be an advantage 
in using the waveform parameters, rather than the source 
parameters that naturally come out of rigid body calcula¬ 
tions. Let us look at the waveform parameter description of 
the gravitational wave signal for the three particular cases 
of interest. We have two goals: (i) to relate the source pa¬ 
rameters to the waveform parameters, and (ii) to identify 
sensible ranges to search over in the waveform parameters. 
A summary of the identified parameter ranges is given in 
the Appendix; see table 


3.1 Triaxial star, not spinning about a principal 
axis 


Starting with equation (46 1 , the motion of a triaxial non- 


aligned star leads to (Jones 


20121 


t 22 

I = —e 


— 2 i(r 2 t + 0 o) 


- [A 721 (sin^ ijj — cos^ ip cos^ 6 ) — 


A/si sin 6 + i A /21 sin 2 ip cos 6], (58) 


t21 

I = —e 


so that 


— 2 (nt+ 0 o) 


^ [A 721 sin 2 ip sin 6 -|- 
0 


i(A 72 i cos^ Ip — A/ai) sin 20 ], 


(59) 


c: 


c: 


complex _ ^[A72l(sin^ l/) — COS^-i/lCOS^ 0) 

I t) 

—A/31 sin^ 6 + 2 A/21 sin 2 ip cos 6 ], 

/ StT r 


complex 

21 


- [A721 sin 2i/) sin 0 -|- 
i(A721 cos^ p) — A/si) sin 20], 


(60) 


(61) 


from which we see 

C 22 = —2{[A72i(sin^ V'“ cos^'i/icos^ 0)— 

A/si sin^ Og -|- (A721 sin 2'!/) cos 0)^}^^^, (62) 

C 21 = —2{(A72i sin2'!/)sin0)^-I- 

r 

(A721 cos^ ip — A/si)^ sin^ 20}^^^, (63) 

$22 ~ 2 (()gw,o 


$?i 


A 721 sm2p} cos 0 
^21 (sin^ Ip — cos^ ip cos^ 0 ) — A/ai si 

O, 


sin 


if) 

(65) 


If values are given for the quantities {( 522 , $^ 2 , < 521 , $^i}, 
as would be the case in the event of a detec¬ 
tion, the above four equations in the five unknowns 
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{f2^A72i/r, n^A/si/r, 0, 0 gw,o} would then generate a 1- 
parameter family of solutions. Note that, when evaluating 


the inverse tangent functions of equations (641 and (651, care 


must be taken to select the correct root so as to correctly re¬ 


construct the complex mass numbers of equations (601 and 

(6T|. 


Having found the algebraic relationship between the 
source and waveform parameters we can now turn to the 
issue of selecting ranges in the waveform parameters. Care¬ 
ful study of equations (|62[) and (|63[) shows that if one selects 


A 721 and AI 31 according to equation (221 then the corre¬ 


sponding bounds on the amplitude parameters are: 

0 < C22 < -A7niax, 

r 

2 Q^ 

0 < C 2 I < -A7niax. 

r 

For the two phase parameters the default range is: 


As discussed in Section |2.2[ we are always free to insist 


0 < 6 ^ < 7 r/ 2 , as made explicit in equation (34l. With this 
choice, both sin^ 9 and sin 29 will always be non-negative. 
In contrast, A 731 can be either positive or negative, so we 
should treat the Adai > 0 and Adsi < 0 cases separately. 
For A731 > 0 case, we can read-off 



^A73isin^0, $ 22 = 0 , 

5 


A73isin20, $21 = 7r/2. 


(76) 

(77) 


Converting to the parameters C 2 m and <E> 2 m that actually 


( 66 ) 

appear in 

the waveform: 


2 H^ 2 

(67) 

C 22 ~ 

-A731 sin^ 9, 


r 



2 H^ . 


C 21 = 

-A731 sin 29, 


4*01 — t 


^gw,0 ' 


(78) 

(79) 


0 < 4>L < 27r. ( 68 ) 

Together with the ranges given in equations , these 

ranges are sufficiently wide to cover all physically distinct 
stellar configurations, and accommodate all possible ad¬ 
ditional information obtained by non-gravitational wave 
means. 

However, from the point of view of carrying out a 
gravitational wave search without such extra information, 
these ranges are redundantly large. The polarisation com¬ 
ponents change sign under the operation 4*^ —>■ 4*^ -|- tt. 
The waveform h{t) also changes sign under the operation 
V’poi V'poi + Tr/S, so we can halve the range in one or other 
of the polarisation angle or the phases. We therefore have 
the options: 

0 < V’poi < 7 r/ 2 , 0 < ^ 2 m < Stt, or (69) 

0 < V’pol <71, 0 < ^ 2 m < TT- (70) 

If one or other of these restricted ranges are employed in a 
search, and a detection is made with parameters (V'poi, 4 > 2 m)j 
three other equally acceptable solutions can be obtained 
through successive applications of the transformation 

T^pol —>■ T^’pol + 7r/2, ^ 2 m ^ 2 m + TT. (71) 

See the first column of Table for a summary of these 
possible choices in parameter ranges. 


These equations can be inverted to give: 


n^Ahi 


1 + 

^ C 21 \ 

(80) 

r 

V 2(722 / 

(^gw,0 — 

4*^1 + 

TT 

2 ’ 


(81) 

tan 9 = 

26*22 



(82) 


C 21 


Note that in this case 

4*22 “ 24*21 +7r, (83) 


a relation that could be hardwired into any search for oblate 
biaxial stars using the waveform parameterisation. 



For A731 < 0 case, we can 

read-off 


C 22 

/StT . ■ 2 n 

= -W—A 73 ism 9, 

V 5 

4*22 = tt. 

(84) 

C 21 

= —-\/^A 73 i sin 20 , 

V 5 

4*21 = -n/ 2 . 

(85) 

Converting to the parameters 
appear in the waveform: 

C 2 m and 4 * 2 m that actually 

6*22 

= - —A73isin"0, 

r 

4*22 — 2(()gw^o + TT, 

( 86 ) 

6*21 

20 ^ 

= - — A /31 sin 26», 

r 

^21 = '(’giv.O + 

(87) 


3.2 Intermediate case: A biaxial star, not 
spinning about a principal axis 


Setting A 721 


0 in equations (58 1 and (591 leads to 


r22 -2i{nt+4>n) /StT . , ■ 2 r. 

I = e ^ A 73 ism 9, 

V 5 

7-21 . — i(dt-\-Sn') / 87r . j ^ 

7 = le ^ — A731 sm20, 

V 5 


so that 


y^complex 

^22 “ 

1 /^A/31 sin^ e, 

V 5 

y^complex 

^21 “ 

A/31 sin 26 e'^^^ 

V 5 


(72) 

(73) 


(74) 

(75) 


These equations can be inverted to give: 


n^Ahi 

r 


0gw,O 




1 -b 


4*?i - 


2 ’ 



with 9 given by equation (82 I. Note that in this case 


4* 


C 

22 


24* 


C 

215 


( 88 ) 

(89) 

(90) 


a relation that could be hardwired into any search for pro¬ 
late biaxial stars using the waveform parameterisation. In 
the event of a successful detection, the measured values for 
the waveform parameters could then be inserted into the 
equations above, to deduce the corresponding source param¬ 
eters. 

To identify the ranges in these waveform parameters to 
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search over, we can convert the ranges in the source param¬ 


eters of equations (331 and (341 using equations (78 1 -( 791 


above (or, equivalently, equations (86l-(87l), to give: 


2 H^ 

(91) 

0 < C2m < -A/niax, 

r 

0 < < 2 vr, 

(92) 


being mindful to use both equations ( |83[ ) (for oblate stars) 
and equation (901 (for prolate stars) to calculate $22 as a 


function of 4 > 2 i- Together with the parameter ranges of equa¬ 
tions these ranges are wide enough to accommodate 

all physically distinct stellar configurations. 

However, from the point of view of carrying out a grav¬ 
itational wave search, these ranges are redundantly large. 
The waveform changes sign under the operation -i/ipoi —>■ 
ippoi + vr/2. It also changes sign under the operation $^1 —>■ 
<E>^i + 7r and simultaneously, swapping the relationship be¬ 
tween $21 and < 1?^2 from equation (831 to (901, or vice versa, 
which simply amount to the transformation $22 —>■ ‘I ’22 + tt. 
It follows that, in the absence of additional non-gravitational 
wave information, we can carry out a gravitational wave 
search over the reduced ranges of either 

0 < i’poi < Tr/S, 0 < 4>21 < 27r, (93) 

allowing for both the oblate and prolate relations of equa¬ 


tions (831 (for oblate stars) and (|90|l (for prolate stars) in 


calculating $§ 2 (^ 21 ) 

0 < i/’poi < TT, 0 < <I> 2 i < 27r, (94) 

with only one or other (but not both) of the oblate and pro¬ 


late relations of equations (831 (for oblate stars) and (901 in 
calculating 4 >^ 2 (^’m)- K either of these restricted parameter 
spaces is used in a search, and a detection is made, with 
parameters ('i/'poi, 4 ? 2 i, “I’m), three other equally acceptable 
solutions can be obtained through successive applications of 
the transformation 


V’poi ->■ V’pol + 7r/2, 




^ 2 m + '?'■ • 


(95) 


Alternatively, i n a s earch for a general triaxlal body (as 
described in Section 3.1 1 , finding a relationship between <E>^i 


and $22 of the form of either of equations (831 or (901 would 


be a sign that the detected signal is coming from a biaxial 
star. 

See the second column of Table |C2| for a summary of 
these possible choices in parameter ranges. 

An alternative choice would have been to instead use 


equations (76 1 and (771 for both the A /31 > 0 and A /31 < 0 


cases, with the understanding that C 22 and C 21 can now 
be either positive or negative, but both of the same sign 


(i.e. both positive, or both negative). Equations (78l-(83l 


then apply in both the oblate and prolate cases. The wave¬ 
form can the made made to change sign under the operation 
C 2 m —> —C 2 m (applied simultaneously to both the C 22 and 
C 21 ). It follows that, in the absence of other information, 
one can search over the reduced parameter ranges of either 


2 fl'^ 

0 < 'i/’poi < vr/2,-A/„ 

r 


< C 2 m < 


20 .^ 


AI^ 


(96) 


0 < 'tppoi < vr. 


„ r 

0 < C2m < - Ain 

r 


(97) 
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If using such a reduced parameter ranges, other equally ac¬ 
ceptable solutions can be generated through successive uses 
of the transformation 

v/pol —t v/pol -(- 7r/2, C2m —^ — C2m- (98) 

We mention this possibility as, while not fitting into the 
scheme of C 2 m being the modulus of a complex number, 
It has the advantage of possibly being easier to implement, 
as there is only one relation connecting the phases 4>22 and 
'h^i, regardless of whether the body is oblate or prolate, so 
some users may find it easier to integrate into their search 
method. 


3.3 Simplest case: A triaxial star, spinning about 
a principal axis 


In this case we can set 0 = 0 in equations ( |58| l and ( |59[ ), so 
that = 0 , while 


V 5 


so that 


(99) 


C2cr'’ie" = t/^A/2ie-2'’f. (100) 

V ^ 

If we follow the convention used in Section 12.31 and choose 
to insist that A /21 > 0 , we can then immediately read-off 



‘I'22 = —2ip. 


( 101 ) 

( 102 ) 


Using the relation 
inverted to give: 


of equation (551 these equations can be 


n^A /21 

r 

2 ((/)gw,o + ' 0 ) 


2 ^ 22 , 

^ 22 - 


(103) 

(104) 


Note that the parameters (/)gw,o and ip are degenerate, as 
expected for this case. In the event of a detection, the mea¬ 
sured values of ( 022 ,^ 22 ) could be inserted into the above 
equations to compute the corresponding source parameters, 
which in this case are related In a very straight-forward way. 

We have already identified ranges in the source param¬ 
eters that cover all possible physicality distinct stellar con¬ 


figurations; see equation (411. These immediately translate 


into the waveform parameter ranges: 


on 2 


These parameter ranges, together with those of equations 
(j^-([^, will be wide enough to accommodate all possible 
non-gravitational wave priors. 

However, from the point of view of carrying out a grav¬ 
itational wave search without such prior information, these 
ranges are redundantly large. The waveform changes sign 
under the operation ippo\ —>■ ippoi + vr/2 and also under the 
operation < 1 >^ —>■ $^2 + vr, and so we can reduce the ranges 
in one or other (but not both) of those parameters: 

0 < ippoi < vr/2, 0 < ^22 < 27r,or (106) 

0 < ippoi < vr, 0<^22< vr. (107) 
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The first choice is the one that reflects the choice tradition¬ 


ally made in gravitational wave searches (see e.g. Aasi et al. 
( 2014[ )). In the event of a successful detection with parame¬ 
ters ('i/’poi, '^^2), three other equally acceptable solutions can 
be generated by successive applications the of the transfor¬ 
mation 


V'poi ->■ 'i/’pol + 7r/2, 


^22 


^22 + ■ 


(108) 


See the final column of TablejC^for a summary of these 
possible choices in parameter ranges. 


4 RELATION BETWEEN THE PRIORS 

As described above, in carrying out a gravitational wave 
search one has a choice as to which set of variables are used, 
the source parameters or the waveform parameters. If one 
is using Bayesian methods to conduct the search, one also 
needs to specify prior information on the range of each pa¬ 
rameter, and supply a function giving one’s initial belief as 
to its probability distribution. For the triaxial star rotating 
about a principal axis, the two sets are essentially the same, 
but for the biaxial star, and the triaxial star not rotating 
about a principal axis, there is a non-trivial conversion to 
be made. 

Unfortunately, it is not obvious what a physically mo¬ 
tivated choice of priors would be, in terms of either set 
of parameters. This is particularly true for the amplitude¬ 
like source parameters A72i,A/3i or the wave parame¬ 
ters C 2 i,C 22 - The choice of priors would be related to the 
strength of the solid crust, the precise mechanism producing 
crustal deformation, and, for the model of Jones ( 2010| , the 
strength and orientation of the superfluid pinning. We will 
therefore content ourselves here with a relatively simple con¬ 
sideration: if we make some simple choice of priors for the 
source parameters, we will evaluate the corresponding pri¬ 
ors for the waveform parameters. This will illustrate the fact 
that a simple choice of, say, a relatively simple rectangular- 
type distribution in one set of parameters, does not corre¬ 
spond to such a simple distribution when expressed in the 
other set. 

We will only consider the biaxial case, as there it is easy 
to carry out calculations analytically. The non-trivial con¬ 
version is between the ‘wobble angle’ 6 and the asymmetry 
A/ai for the source parameters, and the amplitudes C 21 and 
C 22 in for the waveform parameters. The relevant formulae 



C 21 = 


C 22 = A/ai sin^ 0. 

5 


(109) 


( 110 ) 


We will take as a simple example of a set of priors the fol¬ 
lowing: 

(i) 9 drawn by choosing a point randomly and uniformly 
from the upper half of the unit sphere 0 < 9 < -k 12. (The full 
range 0 < 0 < tt is not required, because of the degeneracy 
discussed in Section [2.2[ we can always insist our 9 value lies 
in this upper hemisphere). 

(ii) A/ai drawn uniformly over the interval (0, A/max), 
independently of the value of 9, corresponding to an oblate 
star. 


To simplify things, we can work with a dimensionless quan¬ 
tity /: 

A/ai 


I = 


AA 


( 111 ) 


The corresponding separately normalised priors are 

P{9) = sine, 0<e<7r/2, (112) 

P{1) = 1, 0 < 7 < 1, (113) 

giving a joint prior 

P{9,i)=sin9. (114) 

The parameter space is simply a rectangle in (9,1) coordi¬ 
nates, with a probability density that depends only upon 
9. 

To eliminate annoying factors, and made our ampli¬ 
tudes dimensionless, define 


021 = 

021 

(115) 

A7n,ax 

C 22 = 

O 22 

(116) 


^ A r 

5 ^-tmax 

so that our transformation equations become 

C 21 = 7sin2e, 

C 22 = 7 sin^ 9. 


(117) 

(118) 


To see the shape of the parameter space in the 
{ 021 , 022 ) variables we can look at the images of all four 
sides of the rectangle formed by the {9,1) variables: 

(i) The side 0 = 0, 0<7<1 maps to C 21 = C 22 = 0, i.e. 
collapses to the origin. 

(ii) The side 7 = 0, 0 < 0 < 7r/2 also collapses to the 
origin. 

(iii) The side 0 = n/2, 0 < 7 < 1 maps to C 21 = 0, 
C 22 = 7 ^ 0 < C 22 < 1. 

(iv) The side 7 = 1, 0 < 0 < n/2 maps to C 21 = sin 20, 
C 22 = sin^ 0. This can be shown to be equivalent to the 


curve C 21 = +2^(722(1 — 022 )- 

To And the actual probability distribution within this 
closed region we can use the conversion formula 


P(C' 2 i,C' 22 )det(J) = P(0)P(7), 
where J is the Jacobian of the transformation: 

J = 

Then 


/ 0021 
0.1 

dC2i ' 
d9 

\ / sin 20 

72 cos 20 \ 

\ 0022 

dC22 

I \ sin^ 0 

7 sin 20 J 

\ el 

do 

/ ^ 



det(J) 

= 27 sin^ 0 = 

27 ) 22 . 


(119) 


( 120 ) 


( 121 ) 


Eliminating AJai between equations (1091 and (1101 we 
have 


sin 0 = 


1-f 


C 21 \ 
2(722 ) 


2-1 - 1/2 


and so we obtain 

P(( 721 ,( 722 ) = 


[C|i + (2022)^V/^ 


( 122 ) 


(123) 
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Figure 2. Probability density function for the prior probabili¬ 
ties in {C 21 , C 22 ) coordinates corresponding to choice of priors in 
(9, A/ 31 ) described in the text, relevant to a biaxial star. We plot 
the logarithm of the probability distribution function of equation 
||123[l, and truncate the plot close to the (singular) origin. 


A plot showing the prior probability distribution as ex¬ 
pressed in terms of {C 21 , C 22 ) is given in Figure]^ The rel¬ 
atively simple probability distribution of the {8,1) coordi¬ 
nates, non-zero over a rectangular region, has mapped into a 
highly non-uniform probability distribution, bounded by the 
curves described above, with the probability density going 
singular at the origin of the {C 2 i,C 22 ) system. We plot the 
logarithm of the probability distribution to minimise con¬ 
trast between different points, and truncate the plot close 
to the singular origin. This serves to illustrate that a simple 
choice of priors in terms of one set of variables can lead to 
a more complex prior function in terms of the other set. 


5 SUMMARY AND DISCUSSION 


We have examined various models for gravitational wave 
emission from steadily rotating neutron stars. We have 
shown that, when written in term of the conventional pa¬ 
rameters, termed here the ‘source parameters’, the wave¬ 
forms contain a number of discrete degeneracies, with dif¬ 
ferent values of the source parameters giving rise to exactly 
the same received gravitational wave signal h{t). These are 
related to arbitrary choices as to how one lays down the axes 
{Ox, Oy, Oz) on the rotating star. We show how these de¬ 
generacies can be removed by restricting the ranges of the 
source parameters, explicitly giving one (out of a possible 
multitude) of choices for each model. 

We also re-examined the gravitational wave emission 
model of Jones (20101, which described a star containing a 


pinned superfluid. The waveform as written down in [Jones] 
(20101 contained a continuous degeneracy. To remove this 
continuous degeneracy we introduced a new way of writ¬ 
ing down the waveform in terms of so-called waveform pa¬ 


rameters. Any future gravitational wave search for such a 
signal would be best performed using the waveform param- 
eterisation. However, in the event of a detection, it may be 
of interest to convert the waveform parameters back into 
a (1-parameter) family of source parameters. The formulae 
given in this paper will allow such a conversion, and the 
minimal ranges In source parameters identified may help in 
avoiding unnecessary complications in this process. We also 
briefly commented upon the (non-trivial) relation between 
prior probabilities assigned on one set of parameters verses 
priors on the other set. 

The results of this paper will be of use to researchers 
devising methods to detect such mutli-harmoinc continuous 
gravitational wave signals, and in carrying out parameter 
estimation in the event of a detection. A study of precisely 
these issues is underway, and will be presented elsewhere 
(Pitkin, Gill, Jones & Woan, in preparation). 

One other sort of continues gravitational wave search 
that has received attention concerns r-modes. As shown in 


Owen (20101, the gravitational wave signal from a (cur¬ 


rent quadrupole-dominated) r-mode is very similar to that 
of a (mass quadrupole-dominated) steadily spinning triaxial 
star, with a rotation of 7r/4 in ippoi transforming one class 
of signal into the other. It follows that the discussion of tri¬ 
axial aligned stars, of Section [2.3[ applies to r-modes, with 
only minor modification; the quantity A /21 of the triaxial 
aligned star has only to be replaced by a measure the the 
current quadrupole of the r-mode. The size of the current 
quadrupole is often measured by the dimensionless parame¬ 
ter a; see e.g. equation (5) of |Owen (20101. If one mimics the 
convention used here, where we insist A /21 > 0, to require 
a > 0, the r-mode waveform has degeneracies only with re¬ 
spect to the polarisation angle ippoi and the overall phase of 
the signal $gw,o, and the simple degeneracies of section p(3| 
apply again. 

We end with a few comments on electromagnetically- 
derived priors. The spin frequency and sky location, derived 
from radio pulsar observations, are routinely used in tar¬ 
geted gravitational wave searches, see e.g. [Andersson et al(] 
(20111. For a small number of stars, information derived 
from pulsar wind nebulae has been used to provide values 
for the spin inclination angles (t, ippoi). The pulsar wind neb¬ 
ulae observations are in fact insensitive to the difference be¬ 
tween the spin vectors fl and —fl, which in the notation 
used here means they are insensitive to the transformation 
(t, V’poi) ('^r —V’poi+7''), so the value of ippoi is determined 
only up to the addition of multiples of tt. This is sufficient 
to break the discrete tppoi —>■ V’poi + 7r/2 degeneracy con¬ 
sidered in this paper, and so is sufficient break the discrete 
degeneracy in the remaining source parameters. This means 
that a successful gravitational wave detection that made use 
of the electromagnetic information on ijjpoi would provide 
unique values for all parameters (aside from the continuous 
1-parameter degeneracy if the triaxial non-aligned model is 
being considered in terms of source parameters), apart from 
ippoi, which would have the tt ambiguity. 

It is difficult to see how prior information on other pa¬ 
rameters could be obtained, with the possible exception of 
the parameter 4>gw,o. This parameter essentially gives the 
rotational phase of the star’s mass quadrupole at (retarded) 
time t = 0. In principle, one could try and use some the¬ 
oretical model to relate the rotational phase of the star’s 
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mass quadrupole to some electomagnetically-observed fea¬ 
ture. There are (at least) two types of observation that one 
conld potentially use in this regard. Firstly, models where 
magnetic strains are responsible for deforming the star and 
generating the non-zero time varying mass quadrnpole (see 
e.g. Mastrano et al. (20131) might lead to a prediction for the 


relative phasing of the gravitational wave and pulsar signals, 
giving a prediction of the value of <(>gw,o- Secondly, |Bildst^ 
(19981 proposed that temperature asymmetries in the star’s 
crust could generate a non-zero mass quadrupole (see also 


Ushomirsky et al. (2000l). The observation of such a tem¬ 


perature asymmetry (as a modulation in surface brightness) 
could also be used to predict the absolute the phasing of the 
gravitational wave signal, although (as far as we are aware) 
there are no candidate objects of this class of interest for 
gravitational wave searches. 

Note that in both cases it is necessary to ensure the 
model is sufficiently detailed to give the relative phasing 
of gravitational wave and electromagnetic signals, taking 
proper acconnt of the fact that the gravitonal wave signal 
is sourced by an integral over the whole body’s mass dis¬ 
tribution (as per eqnation (46l), while the electromagnetic 
signals will be produced either at or above the surface. In 
the case of making use of the pulsar radio phase, this is 
particularly problematic, as the radio pulsations are likely 
produced at some considerable altitude above the surface, 
and it is not clear if the large scale external magnetic field, 
connected with the pulsations, will be aligned with the large 
scale internal magnetic field, sourcing the mass quadrupole 
deformation. 

Rather than attempting to supply a prior on the ro¬ 
tational phase ()>gw,o, a more pragmatic approach might be 
to be leave this phase unconstrained, and wait until a suffi¬ 
ciently large number of gravitational wave observations have 
been made. A statistical analysis could then be used to see 
if this phase correlates with, say, the pulsar phase. A cor¬ 
relation of this sort would then be a sign (but not a proof) 
that the magnetic fields do indeed play a role in signifi¬ 
cantly deforming the star, potentially providing a discrim¬ 
inant between magnetic-induced deformation and crustal 
strains, braking this degeneracy, and adding significantly to 
the physical insight provided by the gravitational wave ob¬ 
servations. Clearly, more detailed modelling is required to 
properly explore this interesting issue. 
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APPENDIX A: ROTATIONS DESCRIBED IN 
TERMS OF EULER ANGLES 

It will be useful to describe rotations about one of the body 
axes {x, y, z) in terms of transformations of the Euler angles 
{9, (j), Ip). To do so, it is useful to write down the unit vectors 
along the body axes (ej, ey, eg) in terms of their components 
with respect to the inertial frame unit vectors. This is most 
easily done by constructing the matrix that carries out an 
active rotation from inertial axes to body axes. We need to 
first perform a rotation of pj about the inertial 2 -axis, then a 
rotation of 6 about the inertial a;-axis, and finally a rotation 
of (p about the inertial 2 -axis. In an obvious notation, we 
then have 
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Rab = [RjcP)R^{e)Rm;)]ab, 


The relevant active rotation matrices are: 


RlbW 

Rab{9) 

KbW 


cos Ip — sin Ip 0 
sin ip cos ip 0 

0 0 1 

1 0 0 

0 cos 9 — sin 9 

0 sin 9 cos 9 

cos (p — sin 0 0 

sin (p cos <p 0 

0 0 1 


(Al) 


(A2) 

(A3) 

(A4) 
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Allowing Rab to act on the inertial unit vectors (ex,ey,ez), 
we obtain 


es 


cos (f) cos Ip — sin (p cos 6 sin tp 
sin (p cos tp + cos (p cos 6 sin tp 
sin 0 sin ip 

— cos 4 > sin Ip — sin (p cos 6 cos ip 

— sin (p sin ip + cos (p cos 6 cos ip 

sin 0 cos Ip 

sin (p sin 9 

— cos (p sin 9 

cos 9 


(AS) 

(A 6 ) 

(A7) 


We can then see how these body-frame basis vectors trans¬ 
form under operations of the form {9, cp, ip) —>■ {9, cp, ip) in 
order to identify the nature of the corresponding rotation of 
the body axes. 


A1 Rotation of tt about Oz. 


Under 'ip 

Gx 

Ip TV, we see that 

~Gx5 

(A 8 ) 

By 


(A9) 

Bz —^ 

H-Gz- 

(AlO) 


This is clearly a rotation of tt about Oz. 


A2 Rotation of tt about Oy. 

Now consider the transformation 


of 

1 

t 

(All) 

(p —>• (p + TT, 

(A12) 

Ip —>■ —Ip. 

(A13) 

This produces the transformation 

(A14) 

By 

(A15) 

Bz —Bz. 

(A16) 


This is clearly a rotation of tt about Oy. 


the various sine and cosine functions instead: 


sin 9 = 

(l-sin^6»sin^V’)^''^ 

(A23) 

cos 9 = 

sin 9 sin ip 

(A24) 

sin (p = 

— cos (p cos Ip + sin rp cos 9 sin ip 
(1 — sin^ 0 sin^'!/)) i /2 

(A25) 

COSTp = 

sin (p cos Ip + cos (p cos 9 sin ip 
(1 — sin^ 0 sin^ ^)i /2 

(A26) 

sinip = 

cos 9 

(1 — sin^ 9 sin^ ip)''-/^ ’ 

(A27) 

cos Ip = 

sin 9 cos Ip 

(A28) 

(1 — sin^ 0 sin^ ip)^/‘^ 


APPENDIX B: THE TRANSFORM 

.0,20 0,20 

As described in Section it is useful to identify parameter 
transformations such that —>■ as these com¬ 

bined with transformations of the form ipp^i —>■ ippoi + tt(2 
give rise to a to give a symmetry in the waveform, allowing 
for a reduction in the parameter space range. 

As described in Section in terms of the waveform 
parameters, the transformation is simply 

'l‘ 2 m 'l‘ 2 m + • (Bl) 


In terms of the source parameters, in the case of a triax- 
ial aligned star, the transformation can be achieved straight¬ 
forwardly, by the mapping (</>gw,o + V') (^gw,o + V’) + 7''/2, 
as described in equation (45l. In the case of a biaxial star, 
the transformation is again straightforward, as given by the 
mapping A /31 —>■ —AI 31 of equation (38 1 , corresponding to 
a transformation of the form 5p —^ —Sp. 

However, in the case of a triaxial non-aligned star, our 
choice of following the convention common in ridge body 
dynamics of setting I 3 > I 2 > Ii makes the transformation 
5p —>■ —5p more difficult to describe. Suppose the original 
star density perturbation Sp produces perturbations in the 
moment of inertia tensor 5 I\, 5 l 2 ,Sl 3 , so that the moment 
of inertia tensor, referred to the body axes, is 


h 

= /o + SI\, 

(B2) 

h 

= /o + SI 2 , 

(B3) 

h 

= /o -1- SI 3 , 

(B4) 


A3 Rotation of 7 r /2 about Oy. 

For a rotation of 7 r /2 about Oy, we require 


es ^ - 02 , (A17) 

Oy —> TSy, (A18) 

02 —^ -(-Ox. (A19) 

To see what transformation in the Euler angles produces 
this, we need to solve the equations 

esi{9,f,ip) = -es.{9,(p,ip), (A20) 

ey{9,f,ip) = +ey(9,(p,ip), (A21) 

ei{9,f,ip) = +esi{9,(p,ip). (A22) 


where we choose the axes such that I 3 > I 2 > Ii. We are 
always free to make this choice, and it means that our asym¬ 
metry parameters then satisfy A /31 > A /21 > 0, such that 
the largest principal part of the moment of inertia tensor lies 
along Oz, and the smallest along Ox. Given this convention, 
we can’t simply flip the signs of A /31 and A/ 21 . We need 
to find a more complicated transformation, consistent with 
out convention I 3 > I 2 > I\. 

Under the operation Sp —>■ —Sp, the moment of inertia 
tensor transforms to 

h = lo-Sh, (B5) 

1 2 = IQ-SI 2 , (B 6 ) 

1 3 = Iq-SIs. (B7) 


The solution does not correspond to a simple translation or The Ox axis is now the axis of largest principal moment of 

sign change for the Euler angles, so we give the results for inertia, and the Oz the smallest. To put this star into our 
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conventional form we need to rotate by ±7r/2 about the Oy 
axis. We also need to choose the new perturbations in the 
moment of inertia tensor to be 


Sh = 

-Sh, 

(B8) 

Sh = 

-Sh, 

(B 9 ) 

Sh = 

-Sh, 

(BIO) 


so that the new amplitude parameters are now related to 
the old by 


A 721 = 

SI2 - Sh = -Sh + Sh = A/31 - A/21, 

(Bll) 

A/si = 

5/3 — Sli = —SIi -|- Sh = A/31. 

(B 12 ) 

Note that 

with this choice we have A/31 > A/21 

> 0, as 

required by our convention. If we then carry out the trans¬ 
formation 


-t 0 ,^, 4 ’), 

(B 13 ) 

A/21 

—>■ A/21 = A/31 — A/21, 

(B 14 ) 

A/31 

—>■ A/31 = A/31 

(B 15 ) 


A 3 


where the new Euler angles ( 9 , (p, tp) are those of Section 
equations (A 23 l-(A 28 l, appropriate to a 7r/2 rotation about 


Oy, we find h' 


+ ,x 


— h' 


+ ,x > 


as required. 


APPENDIX C: SUMMARY TABLES OF 
PARAMETER RANGES 
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Table Cl. Table giving one possible choice for the ranges in the source parame¬ 
ters (0, ^gw,o> ^-^215 ^^ 3 i) that uniquely label the orientational of any star’s mass 
quadrupole. The ranges in spin frequency, sky location and inclination angle, common 
to all models, are given in equations Note that in the triaxial aligned case, the 

parameters ^gw,o ^-nd i/j are degenerate, appearing only as the sum ((/>gw,o + ^)- the 
case where electromagnetic observations have not provided a value for V’poh the range 
in this angle can be restricted from the (0,7r) interval of equation ^ to the interval 
(0,7r/2). If this is done, then, in the event of a gravitational wave detection, three other 
equally acceptable solutions can be obtained, each differing by the transformations given 
in equations ( |28[ l and for the triaxial non-aligned case, equation ( |38[ l for the 

biaxial case, and equation psf for the triaxial aligned case. 



Triaxial non-aligned 

Biaxial 

Triaxial aligned 

e 

(0,71/2) 

(0,71/2) 

e = o 

^gw,0 

(0, 271) 

(0,27i) 

(0,7l) 

V) 

(0,71) 

Not present 

Degenerate with ^gw,o 

A /21 

(0,A/3i) 

A721 = 0 

(0, A/max) 

A /31 

(A/ 2 I 1 A/niax) 

( A/max} A/max) 

Not present 


Table C2. Table giving one possible choice for the ranges in the waveform parameters 
{C'22? (^2? ^^ 2 ’^ 2 i}’ consistent with any physical stellar configuration. The ranges in 
spin frequency, sky location and inclination angle, common to all models, are given in 
equations 0-0- In the case where electromagnetic observations have not provided 
a value for ippoi, the range in this angle can be restricted from the (0 , tt) interval of 
equation 0 to the interval (0, '7r/2). In the event of a gravitational wave detection, three 
other equally acceptable solutions can be obtained, each differing by the transformations 
"^pol '^pol = '^/2, ^ 2 m ^ 2 m described in Sectionj^ 



Triaxial non-aligned 

Biaxial 

Triaxial aligned 


( 0 , A/max) 

(0, A/max) 

(0, A/max) 


(0, A/max) 

(0, A/max) 

II 

0 


(0,271) 

= ^^21 ^ (oblate stars) 

$^2 = (prolate stars) 

(0,27i) 


(0,27i) 

(0,27i) 

Not present 
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